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OBSERVACIONES AL ALUMNO 

Al comienzo de cada texto encontrará lo siguiente. 

1. UNA PRESENTACION DE NUESTROS OBJETIVOS en la unidad, la cual puede usar para comprobar si ha completado el 
curso satisfactoriamente. 

2. UN DIAGRAMA ESTRUCTURAL que le muestra cómo están relacionadas las diferentes secciones de la unidad. El diagrama 
consta de: 
a) cuadros rojos que muestran la secuencia principal. 
b) cuadros negros punteados (a la izquierda del rojo) que le muestran lo que esperamos ya sepa antes de fijarse un estudio 

detallado. 

e) cuadros negros (a la derecha del rojo) que le ayudan a escoger temas que puede posponer u omitir si tiene poco tiempo. 

3. UN GLOSARIO de las definiciones, notaciones, etc., usadas en la unidad como referencia inmediata. 

4. UNA BIBLIOGRAFIA de los libros que puede consultar cuando haya acabado el trabajo de la unidad. 

B. CONCLUSIONES Y RESUMEN del trabajo realizado al final de cada texto. 


NUMERACION E INDICADORES 

Dentro de cualquier sección nos referimos a los temas por medio de un número. Al referirnos al mismo tema en otra sección, citamos 
cuatro números: número de la unidad, número del capítulo o parte, número de sección y número del tema. Por ejemplo, “Ecuación 3” 
si aparece en la misma sección, y “Ecuación 4.3.1.3”, si está en otra sección. Los tres primeros números se citan en la parte superior 
de la página. Creemos que este es un sistema claro y fácil para consultar cualquier tema. Los indicadores serán una guía valiosa 
y le ahorrarán tiempo. ; 

En la margen derecha del texto hemos tratado de indicar la clase del material. Los términos usados son tan claros que no necesitan 
explicación. Por ejemplo, “Discusión” significa que una vez que haya entendido el tema, no necesita de nuevo recurrir a esta parte 
del texto. Estos INDICADORES tienen por objeto ayudarlo a localizar rápidamente las partes importantes cuando esté repasando 
o tenga que omitir parte del texto por falta de tiempo. 

También hemos tratado de indicar la importancia de los temas por el sistema de estrellas. 

*** Esto es muy importante. Usted debe entenderlo claramente. En algunos casos lo usará con tanta frecuencia, que se familiari- 
zará con ello sin necesidad de memorizarlo deliberadamente. El punto principal de este curso es no aprenderse los datos de me- 
moria; nosotros creemos que es más importante que entienda el trabajo y sea capaz de aplicarlo en la solución de los problemas. 
** Este es de menor importancia, pero tiene material al cual nos referiremos más adelante. 

* Esto generalmente no lo guía a trabajos adicionales en este curso y lo puede omitir si tiene poco tiempo. 

M Indica el final de cada ejemplo, ejercicio o solución. 
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Objetivos 


Después de trabajar esta unidad debe estar en capacidad de: 


(1) 


(ii) 


(111) 
(iv) 


(v) 
(vi) 


(vii) 


(viii) 
(1x) 


1v 


explicar el significado de los términos siguientes: función com- 
plementaria, coeficientes constantes, ecuaciones del movimien- 
to, vibración forzada, vibración libre, operador diferencial lineal, 
orden de un operador diferencial lineal, frecuencia, período, 
masa, resonancia; 

formular la ecuación diferencial del movimiento de un sistema 
masa-muelle en el cual la masa se mueve sobre una recta; 


determinar si un operador diferencial dado es lineal o no; 
encontrar soluciones particulares de una ecuación diferencial 
lineal con coeficientes constantes; 

comprender el enunciado del teorema establecido en la página 
19 y aplicar el teorema en ejemplos particulares; 

escribir las soluciones generales de ecuaciones diferenciales 
lineales, dadas suficientes soluciones particulares para gene- 
rar el espacio solución; 

determinar los valores de las constantes arbitrarias a partir de 
condiciones iniciales; 

determinar el período de una función dada; 

resolver ecuaciones de la forma 


Lf(t) =kcospt o ksenpt, 


donde 


es una ecuación diferencial lineal que tiene suficientes solucio- 
nes particulares para generar el espacio solución, y cos pt o 
sen pt E al núcleo de L. 
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Interpretación de 
la solución 
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Glosario 


AMPLITUD 


CONSTANTES 
ARBITRARIAS 


COEFICIENTES 
ECUACION DEL 
MOVIMIENTO 


ECUACION 
DIFERENCIAL 
LINEAL 


ELONGACION 


FRECUENCIA 


FUERZA 


FUERZA EXTERNA 


FUNCION 


COMPLEMENTARIA 


MASA 


MOVIMIENTO 


La amplitud de una vibración es la distancia en- 
tre la posición extrema de la partícula vibrante y 
la posición de equilibrio. 


Las constantes arbitrarias son números usados 
para denotar las distintas soluciones de una ecua- 
ción diferencial. 


Por ejemplo, las constantes A y B de la fórmula 


X:t+H—>Acos 1 + Bsent 
de la solución general de 
XxX +x=0. 


Véase operador diferencial lineal. 


La ecuación del movimiento de una partícula es 
la ecuación diferencial que se obtiene de la apli- 
cación de la segunda ley de Newton. 


Una ecuación diferencial lineal es una ecuación 
diferencial de la forma 


Ef g, 
donde L es un operador diferencial lineal. 


La elongación de una partícula es la distancia di- 
rigida desde algún origen fijo hasta la partícula. 


La frecuencia de un movimiento armónico simple 
es el inverso del período. 


Físicamente, la fuerza es un concepto primitivo 
de la mecánica y sólo puede definirse con refe- 
rencia a la experiencia. Matemáticamente, la fuer- 
za se puede definir parcialmente a través del mo- 
do en que aparece en la segunda ley de Newton. 


Una fuerza externa sobre una partícula es un tér- 
mino de la expresión de una fuerza que es inde- 
pendiente de la posición de la partícula. 


La función complementaria de una ecuación di- 
ferencial lineal 


Lf=g 
es un elemento del núcleo de L. 


La masa de una partícula es la razón entre la 
fuerza ejercida y la aceleración de la partícula. 


Un movimiento armónico simple es el movimiento 


ARMONICO SIMPLE de una partícula cuya ecuación de movimiento es 


vi 


mX"+sxX=0 


26 


26 


27 


31 


10 


OPERADOR 
DIFERENCIAL 
LINEAL 


ORDEN 
PARTICULA 


PERIODO 


POSICION DE 
EQUILIBRIO 


RESONANCIA 


SEGUNDA LEY 
DE NEWTON 


VIBRACION 


VIBRACION 
FORZADA 


VIBRACION LIBRE 


Un operador diferencial lineal es un operador de 
la forma 


L:f—> kof + Y K Df 


cuyo dominio es un conjunto adecuado de funcio- 
nes reales, donde los coeficientes Ro, ki, ..., Rn 
son funciones reales. 


El orden de un operador diferencial lineal es el or- 
den de la más alta función derivada que aparezca 
en ella y cuyo coeficiente no sea cero. 


Una partícula es un cuerpo que puede ser repre- 
sentado aproximadamente como un punto mate- 
mático móvil. 


Una función f de dominio R tiene período h si 


Fx +h)= fx) 


para todo x € R, donde |h| es el menor número que 
tiene esta propiedad. 


La posición de equilibrio de una partícula es la 
posición en la cual la partícula podría permanecer 
en reposo si no la perturbáramos. 


La resonancia es la tendencia de las vibraciones 
a alcanzar amplitudes muy grandes cuando la 
fuerza externa oscila con el período de las vibra- 
ciones libres del sistema. 


La segunda ley de Newton establece que la ace- 
leración de una partícula es proporcional a la 
fuerza que actúa sobre ella. 


Una vibración (u oscilación) de una partícula es 
un movimiento reiterado de vaivén. 


Una vibración forzada de una partícula es una 
vibración tal que la expresión de la fuerza inclu- 
ye una fuerza externa. 


Una vibración libre (o natural) de una partícula 
es una vibración tal que la fuerza está completa- 
mente determinada por la posición de la partícula. 
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Página 
18 


18 


25 


Notación 


Los símbolos se presentan en el mismo orden en que aparecen en el 


texto. 
R¿ El conjunto de los números reales positivos y el cero. 
La masa de una partícula. 
S La tensión de un muelle. 
O La función cero: 
0: x-—>0 
cuyo dominio es un subconjunto de R. 
D El operador derivación. 
s 
o pl 
m 
L Un operador. 
k; Una función real. 
A Una trasformación lineal (morfismo) de un espacio vectorial 
en otro. 
Xx Un elemento de un espacio vectorial. 
0 El elemento cero de un espacio vectorial. 
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31.1 EL FENOMENO DE LA RESONANCIA 
31.1.0 Introducción 


En esta unidad reuniremos diversos hilos de argumentos que han ido 
apareciendo en distintos pasos de este curso y los aplicaremos a la 
discusión de un fenómeno físico familiar: la resonancia. La palabra 
resonancia se refiere, en primer lugar, al sonido: describe, por ejem- 
plo, la propiedad acústica de los cuartos de baño donde tantas perso- 
nas despliegan sus dotes de cantantes. Sin embargo, usaremos la 
palabra en un sentido más general para referirnos no sólo al sonido 
sino también a cualquier clase de vibración, como las vibraciones de 
un motor eléctrico o las de un automóvil que transita por un camino 
en malas condiciones, o las vibraciones eléctricas producidas en un 
televisor o en un radio por las ondas electromagnéticas del trasmisor. 
Describiremos suficientemente las leyes físicas importantes en nues- 
tros problemas para que el texto sea inteligible sin necesidad de co- 
nocimientos previos de mecánica. 


Usted puede comprobar por sí mismo, con gran facilidad, lo que es la 
resonancia. Cuelgue un objeto pequeño cualquiera, por ejemplo una 
llave, del extremo de una cuerda o una cadena de uno o dos pies de 
longitud. Sostenga en su mano el otro extremo de la cuerda y muévala 
regularmente hacia adelante y hacia atrás sobre una línea horizontal 
de unas cuantas pulgadas de largo. Si hace esto lentamente, el objeto 
seguirá el movimiento de su mano; si lo hace rápidamente, el objeto 
se moverá con dificultad; pero, si lo hace a la frecuencia correcta, el 
objeto oscilará con una amplitud mayor que la de su mano, y usted 
se dará cuenta que basta un movimiento casi imperceptible de la ma- 
no para conseguir que el objeto oscile violentamente. 


Un experimento de otra índole pero que ilustra esto mismo se puede 
realizar en el cuarto de baño. Moviendo su mano (o, si se siente con 
suficientes energías, con todo el cuerpo) hacia adelante y hacia atrás 
en el agua de la tina con la frecuencia correcta, podrá conseguir un 
movimiento oscilatorio lo suficientemente grande para derramar toda 
el agua en el suelo. 
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Otro ejemplo consiste en aplicar una pequeña fuerza intermitente, en 
los momentos precisos, al columpio en que se sienta un niño. Conse- 
guirá, así, mecerlo con una gran amplitud. 


En cada uno de estos ejemplos encontramos que un sistema vibrante 
responde fuertemente a una fuerza aplicada desde afuera si esta fuer- 
za varía con el tiempo de una manera correcta. Esta fuerte respuesta 
se llama resonancia. En esta unidad estudiaremos la resonancia ma- 
temáticamente. 


La investigación matemática de la resonancia usará varias de las 
ideas que ya hemos discutido anteriormente en este curso. Sin em- 
bargo, antes de aplicar la matemática en cualquier forma, tenemos 
que construir una representación matemática, o modelo, de la situa- 
ción física que estudiamos. Este modelo matemático toma la forma de 
una ecuación diferencial pero, como esta ecuación es de segundo or- 
den, las técnicas que discutimos en la Unidad 24, Ecuaciones dife- 
renciales I, no se pueden aplicar directamente para obtener una solu- 
ción exacta (o una fórmula). De hecho, como no existe un método 
deductivo que produzca soluciones exactas de todas las ecuaciones 
diferenciales, nos enfrentamos a una situación de “resolver un pro- 
blema” del tipo discutido por G. Polya en How to Solve It*; es nece- 
sario poner en orden toda la información que parezca de interés y 
utilizarla para construir una solución del nuevo tipo de ecuación di- 
ferencial. Por último, tendremos que examinar las soluciones que 
hemos encontrado e interpretarlas para obtener una explicación del 
fenómeno físico de la resonancia. 


*G. Polya, How to Solve lt, Open University ed. (Doubleday Anchor Books 1970). 
Este libro es el libro base del Curso básico de matemática; nos referiremos a él lla- 
mándolo, sencillamente, Polya. 
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31.1.1 Construcción de un modelo de resonancia 


El primer paso en la construcción de un modelo matemático de la re- 
sonancia consiste en la simplificación de la situación física; después 
podremos ver más fácilmente cómo se describe esta situación en len- 
guaje matemático. Los puntos esenciales de la situación física son: 


(i) un sistema que puede oscilar después de una perturbación ini- 
cial, aun cuando dicha perturbación cese (en los experimentos 
sugeridos anteriormente, tenemos el objeto suspendido de una 
cuerda o el agua de la tina), 


(ii) una perturbación externa del sistema (en los experimentos suge- 
ridos, esta perturbación es el movimiento de la mano). 


Z 


<= 
S 
= 
> 


Uno de los sistemas físicos más sencillos de este tipo es el que se 
ilustra en la figura. Consta de un cuerpo que cuelga de un soporte. Si 
el soporte se mantiene fijo, el cuerpo puede oscilar hacia arriba y ha- 
cia abajo, con lo cual satisface la condición (i), y si movemos el cuer- 
po alternadamente hacia arriba y hacia abajo (ya sea directamente 
con la mano o moviendo el soporte) podemos aplicar al sistema una 
fuerza perturbadora oscilatoria, de acuerdo con la condición (ii). El 
sistema que vemos en la figura no es tan artificial como parece. El 
soporte puede representar uno de los ejes de un automóvil, el muelle 
uno de los muelles del vehículo, y el cuerpo la parte del automóvil que 
está soportada por ese muelle. (Esto es una verdadera simplificación 
porque, por supuesto, no podemos dividir el automóvil en cuatro par- 
tes que se mueven independientemente pero nos sirve para ilustrar 
el modo en que podríamos iniciar un estudio de las vibraciones de un 
automóvil.) 


Se puede considerar también un sistema horizontal como el de la figu- 
ra. El análisis de ambos sistemas es muy semejante y conduce a la 
misma ecuación diferencial. 
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31.1.2 Descripción del movimiento 


En la construcción de una representación o modelo matemáticos de 
este sistema físico nos interesa, en particular, el movimiento de un 
cuerpo en dirección vertical. El modelo matemático consta de dos 
partes: (i) una descripción del movimiento, y (ii) alguna ecuación o 
algunas ecuaciones que representan la ley física que determina el 
movimiento. Así, primero consideraremos la descripción del movi- 
miento; esto es, consideraremos cómo la posición del cuerpo depende 
del tiempo. Este es un problema que radica principalmente en la se- 
lección de una notación conveniente. Recomendamos al lector que 
lea los artículos “Notation” y “Setting up equations” en las páginas 
134-141 y 174-177, respectivamente, de Polya. 


Para simplificar la situación tanto como sea posible, supondremos que 
el cuerpo se puede reducir a una pequeña región del espacio, de ma- 
nera que se pueda representar aproximadamente como un punto ma- 
temático. Un objeto así representado recibe el nombre de partícula. 
Como estamos estudiando el movimiento vertical, supondremos que 
la partícula se mueve a lo largo de una recta vertical. Podemos des- 
cribir la posición de la partícula dando su distancia dirigida desde 
algún punto fijo de la recta (el origen de coordenadas). Esta distan- 
cia dirigida se llama elongación de la partícula; la denotaremos por 
la variable x. 


Movimiento 


Por distancia dirigida entendemos que la distancia hacia un lado del 
origen toma el signo positivo y, hacia el otro lado, tome el signo nega- 
tivo. No importa cuál de los dos sentidos se tome positivo; lo esencial 
es establecer los dos sentidos y ser consistentes con la selección he- 
cha. De una manera completamente arbitraria diremos que son posi- 
tivas las distancias hacia arriba (esto es, los puntos que están más 
arriba que el origen tienen valores positivos de x). Del mismo modo, 
tampoco tiene importancia cuál sea el punto que tomemos como ori- 
gen. No obstante, en los problemas de oscilaciones es conveniente que 
el origen sea la posición de equilibrio de la partícula; esto es, la posi- 
ción en la cual la partícula podría permanecer en reposo si no la per- 
turbáramos. 


a Partícula 
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31,1,2 


Discusión 
* $ 


Definición 1 


$* $ 


Definición 2 


* $ 


En principio, podemos describir el movimiento dando la posición de 
la partícula en todo instante. Así, en principio, el movimiento queda 
descrito por una lista de pares ordenados de la forma 


(tiempo, posición); 
en otras palabras, por una función real 
X: tiempo —> elongación. 


La gráfica de esta función se llama gráfica de elongación-tiempo. El 
diagrama siguiente muestra una posible gráfica de elongación-tiempo 
de un sistema oscilante. 


Elongación 


Gráfica de X 


Para medir el tiempo, escogemos una vez más un origen, esto es, un 
instante particular, a partir del cual se medirá el tiempo. El origen 
más conveniente es, con frecuencia, el instante en que se inicia el 
movimiento: con esta escogencia podemos definir cualquier instante 
del desarrollo del movimiento dando el tiempo (por ejemplo, en segun- 
dos) que ha trascurrido desde que comenzó el movimiento. El tiempo 
trascurrido se suele denotar por la variable t. De este modo, el movi- 
miento está descrito por una función real 


Xitx  (teR¿j) 


escogida de manera que la elongación x en el tiempo t es X(t), es 
decir, 


x = X(0). 


Hemos seguido la práctica normal de tomar como dominio de esta 
función el conjunto Ry* aunque, estrictamente, sería más correcto 
usar el intervalo [0, T], donde T es la duración total de las condi- 
ciones que determinan el tipo particular de movimiento que estamos 
estudiando. Escogiendo Ry* como dominio tomamos una situación 
ideal simplificada porque suponemos que estas condiciones persisten 
por siempre. 


Z 


<= 
pa 


Posición 
verdadera 


, 
xt0) 


Posición de —? 
equilibrio £._....... : 
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Ejercicio 1 = 

Exprese la velocidad y la aceleración (razón de cambio de la veloci- 
dad) de la partícula en el tiempo t en términos de las funciones deri- 
vadas de X. | 


Ejercicio 2 


¿Qué descripción, análoga a la del texto, sería adecuada para el movi- 
miento de una pieza giratoria de una maquinaria, como el rotor de un 
torno? LL] 


31.1.3 Aplicación de las leyes de la mecánica 


Como ya hemos decidido describir el movimiento usando la función 
X, el paso siguiente es encontrar toda la información posible acerca 
de esta función. La matemática sola no nos dice nada acerca de la 
función; debemos recurrir a nuestros conocimientos no matemáticos 
sobre el sistema en consideración. 


Aquí, este conocimiento consiste en las leyes de la mecánica que fue- 
ron descubiertas experimentalmente. Una de las principales leyes de 
la mecánica es la segunda ley de Newton, que establece que la acele- 


ración de una partícula cualquiera dada es proporcional a la fuerza 
total que actúa sobre ella. (No entraremos aquí en el problema de de- 


finir fuerza.) En términos de nuestra descripción del movimiento, la 
segunda ley de Newton se puede escribir 


X"(t) = kF(0), 


donde k es la constante de proporcionalidad y F(t) es la fuerza total 
que actúa sobre la partícula en el tiempo t. Cuanto más pesada es la 
partícula, tanto menor es la aceleración que le imprime una fuerza da- 
da y, por consiguiente, tanto menor es k; consecuentemente, el inver- 


1 z ; , 
so z es mayor cuando se trata de partículas más pesadas. Este inverso 


se denota por m y se llama masa de la partícula. Su valor numérico 
para una partícula dada depende de las unidades de medida que se 
utilicen. 


Podemos expresar esa ecuación diferencial en la forma 
mX"(t) = F(t). 


Esta ecuación es la forma usual de la segunda ley de Newton del mo- 
vimiénto pero nada nos dice todavía acerca del movimiento; para eso 
necesitamos saber lo que significa la fuerza total F(t) sobre la par- 
tícula o, al menos, cómo calcularla. La “fuerza total” a que se refiere 
la ley de Newton significa la suma de las fuerzas aisladas que actúan 
sobre la partícula y, en el caso presente, incluye la fuerza de la grave- 
dad y la fuerza ejercida por el muelle. Para comenzar a analizar el 
comportamiento de este sistema supondremos que estas son las úni- 
cas fuerzas que están actuando sobre la partícula. Esto implica, en 
particular, que despreciamos la fricción y también (por ahora) que 
negaremos la posibilidad de fuerzas externas como la de alguien que 
empuje el peso. El movimiento de un sistema vibratorio que no está 
afectado por ninguna fuerza (excepto la de la gravedad) que venga del 
exterior del sistema, se llama vibración libre para distinguirlo del 
movimiento donde actúa una fuerza externa (de alguna clase) y que 
llamamos vibración forzada- 
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Ejercicio 1 
(2 minutos) 


Ejercicio 2 
(4 minutos) 


31.1.3 


Tema principal 
* * 


Definición 1 
* * 


Definición 2 
* * 


Definición 3 
* + 
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Los experimentos muestran que, con una buena aproximación, en un 
sistema pequeño como el que estamos considerando, la fuerza de gra- 
vedad sobre una partícula es constante y la fuerza del muelle varía 
linealmente con su longitud (ley de Hooke). 


Fuerza contractiva 
(hacia arriba) 


Longitud 
del muelle 


Longitud 
natural 


Fuerza ejercida | por el muelle 


Fuerza expansiva 
(hacia abajo) 


Ley de Hooke 


La pendiente de la recta se llama tensión del muelle y la denotaremos 
por s. Á veces, por conveniencia y sencillez, se toma s = 1 pero aquí 
no lo haremos. 


La fuerza total es la suma de la fuerza del muelle y la fuerza hacia 
abajo debida a la gravedad, que es independiente de la longitud del 
muelle. Por consiguiente, la suma también varía linealmente con la 
longitud del muelle, y la pendiente de la nueva recta también es s. 


Fuerza neta hacia arriba 
F(t) 


e 
[. 
e” 


Fuerza total 


Longitud 
del muelle 


Fuerza gravitacional 


La longitud de equilibrio del muelle es aquella en la cual se alcanza el 
equilibrio de las dos fuerzas que actúan sobre el muelle. Denotamos 
esta longitud por lo y la longitud real la denotamos por l. Podemos 
expresar en forma de ecuación la relación indicada en el diagrama 
anterior entre la fuerza neta hacia arriba F(t) y la longitud del muelle: 


F(t) = (1 — lo)s. 


F(t) 


(continúa en la página 9) 


Solución 31.1.2.1 


La velocidad es la razón de cambio de la posición y su valor en el 
tiempo t viene, por consiguiente, dado por la derivada de X en E 


(velocidad en el tiempo t) = X '(t) 


La aceleración es la razón de cambio de la velocidad y su valor en el 
tiempo t viene, por consiguiente, dado por la segunda derivada de X 
en £, 


(aceleración en el tiempo t) = X" (t) | 


Solución 31.1.2.2 


Aquí, en vez de describir una posición sobre una recta, queremos des- 
cribir la orientación del rotor. En lugar de la escala lineal del ejemplo 
del muelle, tenemos una escala que mide ángulos. 


1,5 radián 


FO 


5 radián 


Motor 


0.5 radián 
l 


El movimiento se puede describir mediante una función real 
f:t—— (lectura del indicador) 


con la lectura del indicador en radianes. Esta representación, sin em- 
bargo, tiene la desventaja de que f(t) no solamente salta de 2r a 0 
cuando el indicador pasa (en sentido contrario al del movimiento de 
las agujas de un reloj) por cero, sino que no registra la magnitud (ne- 
ta) de la rotación. Para evitar estas dificultades se suele trabajar con 
la función 


f:t— (lectura del indicador + 2rn) 


donde n es el número neto de veces que el indicador ha pasado por 
cero, supuesto el sentido positivo como el contrario al del movimiento 
de las agujas de un reloj. Esto es, cada vez que el indicador pase por 
cero en sentido contrario al del movimiento de las agujas de un reloj, 
añadimos 27 a las siguientes lecturas y, por otra parte, cada vez que 
pase por cero en el mismo sentido del movimiento de las agujas de un 
reloj, sustraemos 2r. [5] 
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Solución 31.1.2.1 


Solución 31.1.2.2 


De las definiciones de l, lo y X(t), tenemos que (véase el diagrama) 


lo == l S> X(1). 

Í | 

1 lo 

<= 
NT V 
rr A: 10) 
Por tanto, l — lo = — X(t) y la fórmula de F(t) se puede escribir 

F(t) = — sX(1). 


Hemos obtenido así de la mecánica dos ecuaciones que nos ayudan 
a determinar el movimiento de la partícula en el caso en que las úni- 
cas fuerzas son las del muelle y la gravedad. Estas ecuaciones son 


mX"(t) = F(t) (ley de Newton) 


F(t) = —sX(t) (ley de Hooke) 


Estas son dos ecuaciones simultáneas con dos incógnitas: las funcio- 
nes X y F. El modo usual de trabajar con ecuaciones simultáneas 
consiste en eliminar una de las incógnitas. En este caso estamos más 
interesados en la función elongación X, que describe el movimiento 
de la partícula, que en la función fuerza F. En consecuencia, elimi- 
namos F(t) y obtenemos 


mX"(t) = —sX(t) 
o, equivalente*, 
mX"(t) + sX(t) =0. 


Esta última ecuación se llama ecuación del movimiento del sistema. 
Es una ecuación diferencial que restringe la función X que describe 
el movimiento del sistema; además, esta ecuación selecciona, de to- 
das las funciones concebibles, un conjunto mucho menor de funcio- 
nes posibles. Recuerde que todo esto se basa en muchos supuestos 
físicos, algunos de los cuales son: el soporte es fijo; el peso del muelle 
es despreciable; el muelle obedece exactamente la ley de Hooke; no 
hay fricción ni fuerza externa alguna. 


*«Muchos autores usan la notación x para expresar la elongación X(t), í para la de- 
rivada X'(t) (la velocidad) y Y para la segunda derivada X ”(t) (la aceleración); en es- 
ta notación, la ecuación del movimiento es 


mx + sx =0. 
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(viene de la página 7) 


Definición 4 


$" + $ 
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31.2 MOVIMIENTO ARMONICO SIMPLE 31,2 
31.2.0 Introducción 31.2,0 

Una vez que hemos construido nuestro modelo matemático de las vi- Introducción 
braciones libres y que hemos llegado a la ecuación del movimiento ." 


mX"(t) + sX(t) =0 (teRo), 


podemos prescindir de la mecánica y, por ahora, estudiar esta ecua- 
ción como un problema matemático. Más específicamente, el proble- 
ma consiste en encontrar el conjunto solución de la ecuación del mo- 
vimiento; esto es, el conjunto de todos los posibles movimientos del 
sistema cuerpo-muelle. Estas soluciones son funciones, no imágenes. 
Por esta razón, generalmente es útil escribir la ecuación del movi- 
miento como una relación de funciones en vez de imágenes: . 


mX"” +sX =0. 


Como siempre que nos enfrentamos a un nuevo problema, conviene 
que tengamos en cuenta algunas de las ideas de Polya acerca de la 
solución de problemas. 


31.2.1 Búsqueda de algunas soluciones 31,2,1 
COMPRENSION DEL PROBLEMA Técnica 
$" 


¿Cuál es la incógnita? 


Un conjunto de funciones reales ARE —R 


¿Cuáles son los datos? 


mM, S. 


¿Cuál es la condición? 
El conjunto debe comprender todas las soluciones de 
mX"” +sX =0,. 


TRAZADO DE UN PLAN 


¿Lo ha visto usted antes? 


Si lo ha visto anteriormente, debe encontrar muy fácil el resto de esta 
sección. 


¿Conoce algún problema relacionado con éste? 


En realidad, hemos visto problemas semejantes en la Unidad 24, 
Ecuaciones diferenciales I. Nuestra ecuación actual es también una 
ecuación diferencial pero difiere de las que hemos estudiado allí en 
que es de segundo orden; esto es, la función derivada de mayor orden 
que hay en la ecuación es la función segunda derivada. (En la Uni- 
dad 24 todas las ecuaciones eran de primer orden.) 


Aquí hay un problema semejante y que se ha resuelto antes. ¿Podría 
usted usarlo? 
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En la Unidad 24 resolvimos una ecuación diferencial que es semejan- 
te a la nuestra, excepto en que es de primer orden. La ecuación es 


FAR Z=:0, 
y una solución es 
f:x——> exp (—2x + c) (xeR). 
Aunque este resultado no resuelve la ecuación de segundo orden, 


puede resultar útil en algún momento y, por eso, conviene tenerla pre- 
sente cuando busquemos una solución de la ecuación 


mX” +sX =0. 5 


Por ejemplo, puede ser posible reducir el problema de la solución de 
mX" + sX = O al problema de la solución de alguna ecuación de la 
forma f' + af = O, donde a es un número constante. 


¿Podría resolver una parte del problema? 


Podríamos intentar buscar parte del conjunto solución en vez de bus- 
car el conjunto solución completo. Más específicamente, podríamos 
buscar una sola función X que satisfaciera la condición mX” + sX 
= O. Una vez logrado ese objetivo, estaremos en mejores condicio- 
nes de encontrar la solución general (es decir, el conjunto solución 
completo). 


¿Podría reformular el problema ? 
Hay varios modos de escribir la ecuación; por ejemplo, 
s 
X"+-X=0 
m 


DIX+7X=0 


7 ES 
m 


donde D? significa DoD,i.e. D?:X——= X”. 


Esta última forma es interesante: reformulemos el problema con pala- 
bras. Estamos buscando una función X tal que, cuando se derive dos 


, z s 

veces, sea esencialmente la misma, excepto por el factor — —. O, con 
m 

otras palabras, el operador D? aplicado a la función X produce el 


A An e a S 
mismo efecto que la multiplicación por el número — —. 
m 


Siguiendo esa idea, podemos comenzar por examinar algunas de las 
funciones cuyas funciones derivadas conocemos, para ver si alguna 
de ellas tiene esa propiedad. Este examen nos llevaría a seleccionar 
una función X, derivarla una vez y, después, volverla a derivar. El 
modo más fácil de hacer esto es ampliar la tabla de derivadas que 
dimos en la Unidad 12, Diferenciación I (página 57), como se mues- 
tra a continuación: 
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Í E f" 
AE IA 
xx" x——>mx" 7! x+——>m(m — 1)x"72 
X——>EXp x Xx —> exp x Xx > eXp x 

1 
x=——> In x x= -— AAA SS 
Ñ 3 
xH—> sen x X—>>C0S X x——> —sen x 
LH——+608 x x——> —sen x xH—> —COS X 


(Hemos omitido los dominios.) 


Hay una tabla más en la Unidad 12 pero ésta nos servirá por ahora. 
Hay tres funciones en la tabla que son obviamente posibles, a saber 


exp, sen y cos 
Ninguna de ellas se adapta exactamente porque no hemos producido 
el factor — > pero esto podemos hacerlo con una sencilla modifica- 
ción: 

f E Ss 


xH——>exp (ax)  xH—>aexp(ax)  x+——>a? exp (ax) 


x-—>sen(ax) xH——>acos (ax)  x=——>-—a*sen(ax) 


xH——>CcOs (ax) x—> —asen(ax) x+——> —a? cos (ax) 


(Véase también la Unidad 13, Integración II, página 54, donde aparece 
una tabla de algunas funciones primitivas.) 


e 4 , s s 
En el caso de la función exponencial sería necesario que a? = ——, lo 
m 


cual significa que a es un número complejo. Esto no desecha esta po- 
sibilidad pero comprende algunos conceptos nuevos. Veamos las otras 
dos posibilidades. En cada caso, 
az S 
a* =—- de modo que a= +0, donde w= /—, 
m 


y tenemos unas cuantas posibles soluciones de nuestra ecuación di- 
ferencial: 


X :tH—=>cCOs (at) 


0) ) (teR¿), 


X :tt—>sen(at) 


donde a = +uw y o e 
m 


(Puesto que cos(—at) = cos at, tenemos tres soluciones en vez de 
cuatro.) 


¿Puede comprobar el resultado? 


Como hemos obtenido estas soluciones a través de varios pasos explo- 
ratorios, debemos ahora comprobar que estas funciones verifican la 
ecuación original. Hágalo mediante una sustitución directa en la 
ecuación diferencial. Otra comprobación interesante es notar que es- 
perábamos obtener un movimiento oscilatorio y esto es precisamente 
lo que representan las funciones seno y coseno. 


E2 


Ejercicio 1 
Encuentre dos soluciones de la ecuación 


y - po, 
donde el dominio y el codominio de g es R. a 
Ejercicio 2 
Encuentre dos soluciones de la ecuación 

g” - 28 -3g=0, 
donde el dominio y el codominio de g es R. lp] 
Ejercicio 3 
Hemos encontrado unas cuantas soluciones de nuestra ecuación 

Pa 

m 


de la partícula del muelle. ¿Hay algunas otras? ¿Puede sugerir alguna 
más? ¿Ha conseguido usted todas las soluciones? 


(No damos la solución de este ejercicio. En las páginas siguientes tra- 
taremos estas preguntas pero puede ir pensando acerca de algunos de 
estos puntos antes de continuar la lectura.) |] 


31.2.2 Búsqueda de la solución general 


El análisis realizado en la sección anterior nos abre algún camino ha- 
cia la solución general (es decir, el conjunto solución completo) de la 
ecuación mX” + sX = O que hemos estado examinando, pero aún nos 
queda camino por recorrer. Hemos encontrado solamente unas cuan- 
tas soluciones de la ecuación del movimiento pero, para alcanzar un 
conocimiento pleno de las vibraciones libres del sistema físico que 
estudiamos, necesitamos estar seguros que conocemos todas las so- 
luciones. 


Hay varias maneras de acercarnos a este problema de la generaliza- 
ción de unas cuantas soluciones conocidas, de modo que obtengamos 
la solución general, es decir, el conjunto solución completo. Por ejem- 
plo, podríamos usar las ideas descritas en “Generalization”, Polya, 
página 108. Esto nos daría un camino que nos llevaría a la respuesta 
de nuestro problema. (Aborde el problema usando estas ideas antes de 
seguir adelante.) Sin embargo, hay a mano un método de mayor fuer- 
za, que también está relacionado con uno de los artículos de Polya: 
“Auxiliary Problem”, Polya, página 50. Como establece Polya, no exis- 
te ningún método infalible para descubrir problemas auxiliares ade- 
cuados, aunque él da dos sugerencias específicas. No obstante, hay 
un problema auxiliar que nos ha sido útil desde la Unidad 3, Opera- 
ciones y morfismos, en adelante: encontrar el morfismo. Usted recor- 
dará que un morfismo consiste en una función h de un conjunto en 
otro (o en sí mismo) conjuntamente con una operación binaria +. y una 
operación binaria [J, definidas, respectivamente, en el dominio y en 
el conjunto imagen, que verifican la condición 


h(x o y) = h(x) O hy) 


13 


MB 31.2.1, 31.2.2 


Ejercicio 1 
(3 minutos) 


Ejercicio 2 
(5 minutos) 


Ejercicio 3 
(5 minutos) 


31.2.2 


Tema principal 
de de 


(continúa en la página 14) 


Solución 31.2.1.1 
La ecuación se puede escribir como 
D*g =g, 
de modo que D? tiene el mismo efecto sobre g que la multiplicación 
por 1. Esto sugiere la función exponencial. Si probamos t —— exp(Íat), 


encontramos que requerimos que a? = 1, i.e. a = +1. Por tanto, 
hay dos soluciones 


t—e (teR) y tHe" (teR). 
Usted debe verificar que ambas funciones son soluciones de la ecua- 
ción diferencial g” — g = O. E] 
Solución 31.2.1.2 


Puesto que la función exponencial t-—> exp(at) y sus funciones de- 
rivadas son iguales excepto en los factores a que se introducen en ca- 
da paso, tenemos todavía esperanza de ajustar el factor a para en- 
contrar algunas soluciones. Por tanto, probemos 


g :tH—— exp (at) (te R); 
entonces, 


g"(t) — 2g(1) — 3g(t) 


a? exp (at) — 2a exp (at) — 3 exp (at) 
= (a? — 2a — 3) exp (at). 


Esta expresión se anula para todo t si 


a? -2a-3=0 
es decir, si 
(a- Va+1)=0 
o sea, si 
a=30a=-l. 


Así, pues, las dos soluciones son 
g :t—— exp (31) (te R); 
g:tH——>exp (—1t) (te R). 


Verifique que ambas funciones son soluciones de la ecuación dife- 
rencial. ”m 


(viene de la página 13) 


para todos los elementos x y y del dominio de h. Aplicaremos la idea 
de morfismo en un contexto de espacio vectorial, de manera que nos 
ayude a encontrar la solución de nuestro problema presente. 


Probablemente no encontrará muy natural esta sugerencia a causa 
de las ideas novedosas que encierra; sin embargo, esta es, precisa- 
mente, la clase de sugerencia que es básica en el pensamiento mate- 
mático. Trataremos de explicar el proceso dentro del cual se desarro- 
lla esta idea. Estamos intentando resolver una ecuación; lo primero 
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Solución 31.2.1.1 


Solución 31.2.1.2 


que nos interesa cuando resolvemos una ecuación es conocer explíci- 
tamente el conjunto donde está definida la ecuación. En general, las 
ecuaciones se pueden compendiar en la forma 


f(x) = b, 


donde f es alguna función (ocasionalmente es una aplicación que no 
es función pero, en todo caso, la supondremos sencilla), b es un ele- 
mento conocido del codominio de esta función y x (un conjunto de 
valores de x) se ha de determinar en el dominio de f, que es el conjun- 
to en el cual la ecuación está definida. Si b está en el conjunto imagen 
de f, entonces existe una solución (o soluciones); en caso contrario, 
no hay solución. Aunque cada tipo de ecuación requiere, generalmen- 
te, alguna técnica particular, hay algunas ideas que se aplican a mu- 
chos tipos de ecuaciones. Queremos descubrir y separar esas ideas. 
Para ello, tratemos de determinar cualesquiera propiedades especia- 


les de f y de su dominio. 


¿Hay unas operaciones (o relaciones) en el dominio que nos den una 
estructura conocida? (Entendemos por estructura algo así como gru- 
po o espacio vectorial.) ¿Es f un morfismo respecto a cualquiera de las 
operaciones binarias obvias del dominio? Ya hemos visto en la Uni- 
dad 23 y en la Unidad 26, Algebra lineal Il y IT, respectivamente, có- 
mo nos puede ayudar una estructura de espacio vectorial a la solu- 
ción de ecuaciones. También nos puede ayudar una estructura de 
grupo. 

Para decirlo brevemente, la ecuación diferencial que estamos tratan- 
do de resolver es “precisamente otra ecuación”. Su solución compren- 
de algunas peculiaridades porque es una ecuación diferencial pero 
también comprende algunos principios algebraicos generales. Ya he- 
mos visto algunas de las peculiaridades de las funciones soluciones 
que ya hemos encontrado. Ahora buscaremos algunos de los princi- 
pios generales. 


Las ideas comprendidas disminuyen y acortan las fronteras en la ma- 
temática. No continuaremos aislando una ecuación en el tema donde 
tiene lugar, sino que la observaremos desde un punto de vista más 
general. 


Retornemos al caso presente: estamos interesados particularmente en 
el conjunto de todas las funciones que verifican la ecuación mX"” + 
sX = O. La parte característica de nuestra ecuación diferencial es 
la expresión mX" + sX y, por tanto, debemos estudiar el operador 
definido por esta expresión, esto es, 


L:if——omf"+sf  (feP). 


(El dominio, F, de L, será algún conjunto de funciones. Podríamos es- 
pecificarlo con mayor precisión pero, puesto que ahora no nos esta- 
mos preocupando tanto por el rigor, sería una digresión hacerlo.) Po- 
demos escribir nuestra ecuación diferencial en la forma muy concisa 


LX) =0, 


pero esta notación simplificada solamente es útil si las propiedades 
importantes del operador L también se pueden expresar en una forma 
sencilla. Es aquí donde encaja el morfismo. 


Para satisfacer la definición de morfismo, necesitamos operaciones 
binarias en el dominio y en el conjunto imagen del operador L. Las 
operaciones binarias más importantes que se pueden aplicar a las 
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funciones se discutieron en la Unidad 1, Funciones, y son la adición, 
la multiplicación y la composición. Pensemos por un momento cuál 
de estas Operaciones parece más a propósito para producir un morfis- 
mo cuando se usa en L (es decir, cuál es compatible con £). Las re- 
glas de la derivación, dadas en la Unidad 12, Diferenciación I, nos 
dan una pista; solamente en la adición encontramos un morfismo, 
el que se describe por la regla 


(f+ 8) =S' +g. 
Esto sugiere que se pruebe la adición en el dominio y en el conjunto 
imagen. Encontramos que 

LS + 8) = mf + gY + s(f + 8) 

=(mf” +sf) + (mg” + sg) 

o, en otras palabras, 

Lf + 8) = UN + Lg)  (f,geF), 
de manera que hay un verdadero morfismo 

LF, +)—>(L(F), +) 
asociado al operador L. 


Es posible, desde luego, aplicar el morfismo, de una vez, a la solución 
de nuestra ecuación diferencial pero el trabajo será más fácil si pri- 
mero aplicamos un poco de álgebra. Aunque no hemos especificado 
detalladamente el dominio F, podemos suponer, sin dificultad, que 
F es un espacio vectorial respecto a la adición de funciones y a la 
multiplicación de funciones por números reales. La propiedad de mor- 
fismo de L demostrada anteriormente es una de las dos propiedades 
de morfismo que caracterizan un morfismo entre dos espacios vecto- 
riales, es decir, una trasformación lineal (véase la Unidad 23, Algebra 
lineal II); es, pues, natural investigar si L también verifica la otra 
propiedad de morfismo: 


Laf)=aL(f) (ueR y feF). 


Usando las reglas de la derivación podemos comprobar que, efectiva- 
mente, esta ecuación se verifica puesto que 


Llaf) =m(afY + slaf) = mf” + sf) = af). 
Así, el operador L tiene dos propiedades: 
Uf + g) = LS) + Lg) 
Laf) = aL(f) 


y, por consiguiente, satisface la definición de morfismo de espacios 
vectoriales; esto es, L es una trasformación lineal. 


Como hemos demostrado que nuestra ecuación diferencial se puede 
escribir en la forma 


L(X) = O, 


donde L es una trasformación lineal, podemos usar la teoría de las 
trasformaciones lineales discutida en la Unidad 23, Algebra lineal Il, 
como una ayuda en la solución de esta ecuación. Ya hemos estudiado 
ecuaciones de esta forma en la Unidad 23. Al conjunto solución lo lla- 
mamos núcleo de la trasformación lineal, y vimos que tenía algunas 
propiedades notables: 
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(i) el núcleo es un espacio vectorial (un subespacio de F); 
(ii) dimensión de L(F) = (dimensión de F) — (dimensión del núcleo). 


La segunda propiedad no nos ayuda directamente (no hemos especi- 
ficado F) pero la primera es de un valor incalculable; nos dice que el 
conjunto solución de nuestra ecuación diferencial es un espacio vec- 
torial Expresado en lenguaje sencillo esto significa que, si X y Y son 
dos soluciones cualesquiera de la ecuación diferencial, entonces toda 
combinación lineal de X y Y es también una solución de la ecuación 
diferencial. 


Estamos muy cerca de la solución de nuestro problema si encontra- 
mos el conjunto solución de la ecuación diferencial mX" + sX = O. 
Conocemos algunas “soluciones” de esta ecuación, a saber, las fun- 
ciones 


X :tH——>cO0s (at) 


p (teRg¿), 


X :tt——sen(at) 


donde 


m 


y también sabemos que el conjunto solución es un espacio vectorial. 
¿De qué manera nos pueden ayudar estas soluciones conocidas a es- 
pecificar todo el conjunto solución? Vimos en la Unidad 23, Algebra 
lineal II, que el modo más conveniente de especificar un espacio vec- 
torial consiste en dar una base del espacio vectorial, esto es, un con- 
junto de vectores bases, de manera que todo elemento del espacio 
vectorial se pueda expresar, de manera única, como una combinación 
lineal de los vectores bases. Es natural pensar que, quizá, las solucio- 
nes que hemos encontrado formen una base del espacio vectorial que 
constituye nuestro conjunto solución; en otras palabras, que toda so- 
lución de la ecuación diferencial es una combinación lineal única de 
las soluciones que ya hemos encontrado. Hay un inconveniente: una 
base está formada por elementos linealmente independientes, y sola- 
mente hemos encontrado dos soluciones que son linealmente inde- 
pendientes porque 


cos (— wt) = cos (wt) 


sen(—wt) = —sen(wt). 


Así, pues, modificamos ligeramente nuestra conjetura y sugerimos 
que el conjunto solución de mX” + sX = O es el conjunto de todas 
las funciones de la forma 


X:tH——acoswt + Pfsenot  (teRj), 


E y , S 
donde « y $ son números reales arbitrarios, y w =,/—. 
m 


Ejercicio 1 


Verifique, mediante una sustitución directa en la ecuación diferen- 
cial mX” + sX = O, que toda función de la forma anterior es una 
solución de esta ecuación diferencial. 5 
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Ejercicio 1 
(3 minutos) 


Ejercicio 2 


Hemos aceptado que las dos funciones que estamos usando en nues- 
tra base son linealmente independientes. ¿Puede demostrar que, efec- 
tivamente, son linealmente independientes? (Recuerde que dos vec- 
tores v| y v2 son linealmente independientes si la condición 


av, + fu, =0 


implica que los números reales a y f£ son ambos iguales a 0.) pol 


31.23 Un teorema 


Todavía quedan algunos cabos sueltos que debemos reunir antes de 
que podamos asegurar que hemos encontrado la solución general de 
nuestra ecuación mX” + sX = O. Lo que tenemos que hacer es en- 
contrar, mediante un proceso de luminosas conjeturas o, como dice 
Polya, “por medio de un razonamiento heurístico”, un conjunto de 
soluciones de la ecuación. (Véase “Heuristic reasoning”, Polya, 
página 113.) Puesto que no hay un modo obvio de aumentar el con- 
junto, es natural que conjeturemos que este conjunto es, de hecho, 
el conjunto de todas las soluciones. El razonamiento heurístico es 
bueno en sí mismo pero es un error grave confundir el razonamiento 
heurístico con una demostración rigurosa. Otro modo de plantear es- 
to consiste en considerar que hemos resuelto el “problema de encon- 
trar” pero todavía nos queda el “problema de probar”. 


Puesto que el curso básico está dirigido hacia la exposición de los 
principales conceptos matemáticos y no hacia las demostraciones ri- 
gurosas, no entraremos aquí en la técnica de la búsqueda de tales de- 
mostraciones. Mejor citamos el teorema del que rápidamente se pue- 
de obtener la demostración que buscamos. Tal teorema se aplica a 
una Clase de operadores que incluye el operador específico L: f—= 
mf" + sf que hemos estado considerando. Los operadores de esta cla- 
se se llaman operadores diferenciales lineales; son operadores de la 
forma 


Es PH=3 kk, A DN He EDS + Es RDA SÓ Y 


(fe Fr), 


donde h es un entero positivo que llamamos orden del operador, F, 
es conjunto apropiado de funciones, y ko, ki, ..., kn son funcio- 
nes reales que tienen el mismo dominio que tienen las funciones del 
conjunto F,. (Más precisamente, F, es el conjunto de todas las 
funciones f, cuyo codominio es R y cuyo dominio es algún subconjun- 
to de R, tales que la n-ésima función derivada es continua y tiene el 
mismo dominio que f; también se exige que las funciones ko, ...,kn 
sean continuas.) Especificamos que k, no es la función cero (es de- 
cir, que su conjunto imagen contiene números distintos de cero); de 
no ser así, el término k, X D”f sería cero, lo cual indicaría que ha- 
bíamos usado un valor muy grande de n al escribir la fórmula. Este 
requisito se introduce para evitar tal error al escribir el operador 
f——f', cuyo orden es 1, como f——0 X f” + f', que, aparentemente, 
es de orden 2. 


18 


MB 31.2.2, 31.2.3 


Ejercicio 2 
(4 minutos) 


31.2.3 


Tema principal 
X $ $ 


Definición 1 
kh * * 


Definición 2 


* * 


Como un ejemplo, el operador f——mf"” + sf es un operador diferen- 
cial lineal, con n = 2. 


k,:t——>m 
k, :t 0 (teR¿) 


ko :tH—s 


y F; es el conjunto de funciones reales cuyo dominio es Ry* y que 
tienen, como segundas derivadas, funciones continuas con el mismo 
dominio. En este caso, las funciones kn, ..., ko son todas funcio- 
nes constantes. Puesto que cada k representa una función constante, 
L es un operador diferencial lineal con coeficientes constantes, y las 
funciones k,, ..., ko se pueden tratar como factores numéricos en 
vez de funciones. Señalaremos (sin demostración) las propiedades si- 
guientes: 


(1) F, es un espacio vectorial respecto a las operaciones usuales; 


(2) cualquier operador diferencial lineal L de orden n es una trasfor- 
mación lineal (un morfismo de espacios vectoriales) de F, en 
E). 


(No hay nada esencialmente difícil respecto a las demostraciones de 
estas proposiciones pero son tediosas y poco instructivas.) 


Ejercicio 1 
¿Cuáles de los siguientes operadores satisfacen la definición de un 
operador diferencial lineal? 
NT Df — $? 
UH: DY -f 
(iii) H:f—=>D% — (x—2) 
En todos los casos el dominio es F,. | 


Ejercicio 2 
Escoja uno de los operadores del ejercicio 1 que no sea un operador 


diferencial lineal y demuestre, mediante un contraejemplo, que tam- 
poco es una trasformación lineal. 1ES] 


Ya estamos en condiciones de demostrar un importante teorema. 
TEOREMA 


La dimensión del núcleo de un operador diferencial lineal con coe- 
ficientes constantes es igual al orden del operador. 


Otro modo de formular el teorema es éste: el conjunto solución de una 
ecuación diferencial de orden n de la forma Lf = O, donde L es un 
operador diferencial lineal con coeficientes constantes, es un espacio 
vectorial de dimensión n. El teorema también es válido para una am- 
plia clase de casos donde L no tiene coeficientes constantes, pero no 
necesitamos hacer esta generalización aquí. La demostración del teo- 
rema está más allá del objetivo de este curso pero pensamos incluirlo 
en el curso M201 de matemática lineal. 


Usando este teorema podemos completar ahora nuestra solución de la 
ecuación diferencial mX"+sX = 0. 
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Definición 3 


* $ 


Ejercicio 1 
(3 minutos) 


Ejercicio 2 
(5 minutos) 


Teorema 
$* $ + 


Tema principal 


*os oh 
(continúa en la página 20) 


MB 31.2.2, 31.2.3 
Solución 31.2.2.1 Solución 31.2.2.1 
Si 
X(t) = acos wt + Bsenot, 


entonces la derivación nos da 


Xt) = —awsenowt + Bwcos wt 
y 
X"(t) = —aw? cos wt — Bw?senot 
= —u*X(0, 
de manera que 
mX” +sX =0. Ea 
Solución 31.2.2.2 Solución 31.2.2.2 


Para demostrar que las funciones seno y coseno que estamos usando 
son linealmente independientes, tenemos que demostrar que 


acos + fsen=0=>a= $ =0, 

es decir, 
a.cos wt + Bsenwt = 0 para todo teR¿ 
0 =B:=0) 


Un modo de probar esta proposición consiste en considerar algunos 
valores especiales de t. Si tomamos t = 0, la ecuación que ha de ve- 
rificarse se reduce a a X1+8$8X0=0; estoes, a = 0: si toma- 
mos t = 1r/2w, se reduce a a xO0O+8$8X1=0; ¡.e. 6 = 


| 
Solución 1 Solución 1 
(i) Este operador no es lineal; el término f? lo impide. 
(1i) Este operador es lineal. 
(iii) Este operador no es lineal; el término x ——> 2 lo impide. a 
Solución 2 Solución 2 
(i) Si, por ejemplo, f:x H—— 1, tendremos que 
2H($) = AD — $?) = x——> -2 
pero 
HQS) = DAS) —- QfY = x— —4. 
(i1i) Si f:x——1 tenemos que 
2H(f) = AD — (x—2)) = x—— —4 
pero ' 
HQSf) = D'(2f) — (x-—>2) = x= —2, ES 


(viene de la página 19) 


El conjunto solución de esta ecuación es precisamente el núcleo del 
operador diferencial lineal f——mf"” + sf, que es de segundo orden. 
Por consiguiente, según el teorema, la dimensión de este núcleo es 2; 
el conjunto solución de mX"” + sX = O es un espacio vectorial de 
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dimensión 2. Para especificar este espacio vectorial lo único que ne- 
cesitamos es una base, esto es, un conjunto de 2 soluciones lineal- 
mente independientes. 


Ya sabemos que las dos funciones 
tU ——>COS wt 


y p (te RG) 


tt——senwt 


son soluciones y que son linealmente independientes (ejercicio 31.2.2.2). 
En consecuencia forman una base del espacio solución. El espacio 
solución es el conjunto de todas las combinaciones lineales de los 
vectores bases, esto es, el conjunto de todas las funciones X tales que 


X:tt—acos wt + Bsenwt (teR¿), 


A a . s 
donde a« y £ son números reales arbitrarios, y w = ma 
m 


Ejercicio 3 


Aplique el método que hemos usado en la ecuación mX" +sX = 0 y 
encuentre todas las soluciones de la ecuación diferencial 


D'f - f=0, 
donde el dominio de f es R. 


Esto es, 


(i) determine la dimensión del conjunto solución (un espacio vec- 
torial) usando el teorema; 
(ii) encuentre suficientes soluciones linealmente independientes pa- 
ra formar una base de este espacio vectorial; 
(i1i) escriba una fórmula que dé la solución general en términos de 
esas soluciones linealmente independientes. 16] 


Ejercicio 4 
Encuentre la solución general de 
f"-Sf'+4f=0, 
donde el dominio de f es R. (Véase el ejercicio 31.2.1.2.) Ñ 
Ejercicio 5 


¿Qué está mal cuando aplicamos el método del ejercicio precedente a 
la ecuación 


S"+2f'+$=0, 
donde el dominio de f.es R? ES] 


Para terminar esta sección daremos un ejemplo de cómo se puede ha- 
cer buen uso de algunas manipulaciones formales. 


Ejemplo 1 
Si tratamos de resolver la ecuación 


f"+4f +5f=0 
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Ejercicio 3 
(4 minutos) 


Ejercicio 4 
(4 minutos) 


Ejercicio 5 
(4 minutos) 


Ejemplo 1 


(continúa en la página 23) 


Solución 3 
(i) La ecuación es de la forma Lf = O, donde L es un operador di- 
ferencial lineal de grado 2 con coeficientes constantes. Según el 


teorema, su núcleo es, por consiguiente, un espacio vectorial de 
dimensión 2, y este núcleo es el conjunto solución de la ecua- 


ción. 
(ii) En la solución 31.2.1.1 se demostró que D?f = f tiene como so- 
luciones 
a (te R) 
y 
te" (te R). 


Estas soluciones son linealmente independientes. (Si «ae! + Be” 
= 0 para todo t, entonces tenemos, haciendo t = 0 y después 
t=1l,a+848=0 y ae + Be! =0, y estas dos ecuaciones to- 
madas juntas implican que a = $ =0.) 

(iii) Tomando las dos soluciones linealmente independientes dadas 
en (ii) como vectores bases del espacio vectorial, obtenemos la 
solución general: 


at— e!) + Blt— e”) (te R), 
Le. 
tae + Be”' (te R), 
donde «+ y $8 son números reales arbitrarios. | 
Solución 4 
Si probamos 
Jt—>exp(et) (te R) 
en la ecuación, el miembro de la izquierda se convierte en 
f" —Sf' + 4f = t—(c?e — Sce” + 4e)  (teR) 
= (0? — Sc + 4)(t— e”) (te R). 


Esta ecuación se verifica si el miembro de la derecha se anula; esto 
es, si 


c?—-5c+4=0 
16: 
c=16 4. 
Así, la ecuación diferencial tiene las dos soluciones 
H—=>expt y tH=>exp4£ 


que podemos verificar que son linealmente independientes. Usándo- 
las como vectores bases del espacio vectorial bidimensional de las 
soluciones, obtenemos la solución general: 


t—aexpt + fPexp 4t (te R), 
donde « y 8 son números reales arbitrarios. | 
Solución 3 
Aplicando el mismo método que antes, encontramos que la condición 
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Solución 3 


Solución 4 


Solución 5 


para que ¿t+H——>exp (ct) sea una solución es 
c2+2c0+1=0. 


Esta ecuación cuadrática tiene como única solución c = —1, y por 
consiguiente la ecuación diferencial tiene solamente una solución de 
la forma tt—>exp (ct), a saber, 


t——> exp (—1) (te R). 


Sabemos, sin embargo, que el espacio solución es bidimensional; por 
consiguiente, la función anterior es insuficiente para ser base del es- 
pacio solución. En este, por tanto, el método que hemos usado es ina- 
decuado. Podemos, no obstante, salvar la dificultad; convendría que 
lo intentara. 7 


(viene de la página 21) 


por el método descrito en el texto, encontramos que t -——exp (ct) es 
una solución, supuesto que 


c?+44+5=0 
1.e. 
c=-2+i1 


Esto podría considerarse como un pequeño contratiempo, al menos 
por dos razones: 


(i) siempre hemos supuesto que nuestras funciones eran funciones 
reales; 

(ii) nunca hemos discutido la derivación de funciones que no son 
reales, de manera que, aun cuando aceptáramos que hemos en- 
contrado dos soluciones, no tenemos modo de verificarlas. 


Por consiguiente, podríamos concluir que, si la ecuación tiene solu- 
ciones reales, no son de la forma exponencial. Eso es estrictamente 
correcto, pero no es necesario botar al niño conjuntamente con el agua 
del baño. En la Unidad 29, Números complejos II, llegamos a la fór- 
mula de Euler: 


e? = cos d + ¡send. 


Usando ese resultado podemos escribir 


e 2+M = ¿72 git = e (cost + ¡sent) 


el" 2-M = ¿"UH git = e (cost — ¡sent). 


Ahora, aun cuando la derivación de la función exponencial compleja 
esté en duda, el álgebra sigue siendo válida; el espacio vectorial ge- 
nerado por las funciones exponenciales complejas t-——>exp(—2 + ¡)t, 
tt——exp(—2 — i)t, se puede generar por medio de cualesquiera otros 
dos vectores linealmente independientes que pertenezcan al mismo 
espacio. Nótese que 


e E A O 


2+ il 


e — e. 2-0MM — lie” "sent. 


Olvidándonos de 2 y de 2; (que no son más que multiplicadores esca- 
lares, ya que, en este caso, el conjunto de los escalares es el conjunto 
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ye 


de los números complejos), vemos que 
te cost y tte "sent 


pertenecen al mismo espacio vectorial bidimensional. Se puede de- 
mostrar que son linealmente independientes y que, por tanto, forman 
una base del espacio vectorial. Por consiguiente, usando el “razona- 
miento heurístico”, podemos conjeturar que el conjunto solución de 
la ecuación es 


tae cost + fe" "sent (te R), 


donde « y 6 son números reales arbitrarios. Verifique por sustitu- 
ción que esto es correcto. La independencia lineal significa entonces 
que hemos encontrado el conjunto de todas las soluciones reales de 
la ecuación. E 


Este ejemplo ilustra lo que es un método válido de solución de tales 
ecuaciones. Es interesante notar que estas ecuaciones diferenciales 
no son ficciones de la imaginación de un matemático demente. Tales 
ecuaciones ocurren en la teoría de los circuitos eléctricos, por ejemplo. 
El término e” tiene un efecto de “amortiguación” en el resto de la 
solución; cuando t crece, entonces e” * decrece de modo que, para 
valores “grandes” de t, la solución llega a ser, efectivamente, 


tH=>0 (te R). 


nia 
=! 
ig 
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La función anterior puede aparecer cuando, por ejemplo, se separa un 
sistema de su posición de equilibrio y, después, se deja a sus propios 
recursos, sin ninguna interferencia externa. Esa separación inicial 
producirá oscilaciones que irán muriendo gradualmente, de la misma 
manera que ocurren cuando se lanza una piedra al agua. 


31.2.4 Interpretación de la solución 


En esta sección interpretaremos la solución general 
X :tH——agcos wt + Bsenwt (te R), 


de la ecuación del movimiento de un sistema cuerpo-muelle, en tér- 
minos de los posibles movimientos del sistema. Toda solución de la 
ecuación, donde la amplitud sea menor que l¿, corresponde a un po- 
sible movimiento y, por tanto, la ecuación general que hemos encon- 
trado, con esta restricción sobre a y f, corresponde al más general 
movimiento posible. Podemos usarla para encontrar las propiedades 
generales del movimiento de un sistema vibrante. 
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Discusión 
* * 


31.2.4 


Discusión 
hh 


Lo más importante de estas propiedades es que el movimiento repre- 
sentado por nuestra solución es periódico: se repite exactamente igual 
una y otra vez. Por ejemplo, aquí tenemos las gráficas de las dos so- 
luciones que hemos estado usando como vectores bases de nuestro 
espacio vectorial de soluciones. 


> sent 


” 


Cos wt 


Aquí, T está definida por la ecuación 


wT = 2x1 
Le. 
Ta 
w 


de modo que la primera función, t-——=sen wt, aplica T a sen 27 y 
la segunda, t -— cos wt, aplica T a cos 27. En ambos casos, las 
gráficas de sección comprendida entre 0 y T se repite exactamente 
en las secciones de T a 2T, de 2T a 3T y así sucesivamente. Exacta- 
mente lo mismo sucede en las combinaciones lineales de las dos fun- 
ciones; por ejemplo: 


2senwt-cos wwt 


El número 7 es el período del movimiento. En este caso T depende 
únicamente de la masa del cuerpo y de la tensión del muelle, y no del 
movimiento particular del sistema vibrante. Que el período permanez- 
ca constante es de suma importancia en los sistemas vibratorios. Por 
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Definición 1 
* 


ejemplo, asegura que el tono de la nota emitida por un diapasón, que 
está determinado por el período de las oscilaciones, no se ve afectado 
por la amplitud (el “tamaño”) de las vibraciones, una cualidad muy 


deseable para que el diapasón pueda llenar su función de producir | 


notas patrones. 


Otra característica de la solución general, que podemos relacionar con 
las propiedades del sistema mecánico, es la presencia de dos núme- 
ros no especificados y y fg, que generalmente llamamos constantes 
arbitrarias, en la solución 


X :tH—>acos wt + PBsenot (te R). 


En la Unidad 24, Ecuaciones diferenciales 1, encontramos que la so- 
lución general de las ecuaciones diferenciales de primer orden con- 
tenían solamente una constante arbitraria. Las ecuaciones dife- 
renciales de segundo orden de tipo lineal (que es el único tipo que 
consideramos aquí) tienen una solución general que contiene dos 
constantes arbitrarias. De hecho, una ecuación diferencial lineal de 
orden n tiene, normalmente, n constantes arbitrarias en su solución 
general. 


Para poder escoger de la familia de soluciones (el espacio vectorial) 
la solución particular que describe algún movimiento particular del 
sistema mecánico, necesitamos dos datos informativos que nos per- 
mitan determinar los números « y $. En los problemas de mecánica, 
como el del sistema cuerpo-muelle, esta información generalmente se 
refiere a la posición y velocidad iniciales del cuerpo. Recuerde que 
X(t) es su distancia dirigida, medida desde la posición de equilibrio, 
en un tiempo t, y que X'(t) es su velocidad en el tiempo t. Por ejemplo, 
si el sistema parte del reposo, en el tiempo 0, desde la posición xo, 
entonces tenemos que 


X(0) = xy, puesto que la posición inicial es xp, 


X'(0) = 0, puesto que la velocidad inicial es 0. 


Posición 
verdadera 


Posición de 
equilibrio | A 


Estos dos datos informativos se pueden utilizar para determinar los 
valores de a y $; puesto que, en este caso, X(0) = a y X'(0) = 6, 
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Definición 2 


* * 


obtenemos a = xo, B| = 0. Así, la solución, si el sistema parte del 
reposo desde la posición xo, viene dada por 


XX 9 COS wÉ (teR¿). 


Ejercicio 1 


Si la partícula sale, en el tiempo 0, de la posición xo con velocidad 
Vo, Encuentre la solución apropiada que describe la elongación. Mi 


Ejercicio 2 


Si usted quisiera aumentar la frecuencia (el número de oscilaciones 
por segundo) de un sistema vibrante, entonces, ¿debería 


(a) aumentar la masa? 

(b) poner tirante el muelle? 

(c) disminuir la masa? 

(d) disminuir la tensión del muelle? 


31.3 VIBRACIONES FORZADAS Y 
RESONANCIA 


31.3.1 Un modelo matemático de la resonancia 


Para tratar matemáticamente el fenómeno de la resonancia, debemos 
extender nuestro modelo matemático para que incluya los efectos de 
otras fuerzas distintas de la fuerza —sX (t) que surge directamente de 
la separación del cuerpo de su posición de equilibrio. Las vibraciones 
que ocurren en presencia de tales fuerzas adicionales se llaman vi- 
braciones forzadas, para distinguirlas de las vibraciones libres que 
tienen lugar cuando no están presentes esas fuerzas. Al estudiar la 
resonancia nos interesa particularmente el caso en el cual esa fuerza 
adicional es periódica; ya veremos que las oscilaciones provocadas 
por ese tipo de fuerza pueden alcanzar una amplitud muy grande aun 
cuando la fuerza misma sea pequeña. Un modo de trasmitir una de 
esas fuerzas adicionales al sistema podría hacerse mediante un mo- 
vimiento periódico del soporte. 


Para el tratamiento matemático de esta situación física, supongamos 
que la fuerza adicional tiene la forma F, cos pt, donde F, y p son 
números que podemos usar para ajustar su amplitud y su período. 
También se podría usar un seno en vez del coseno. La razón de usar 
un seno o un coseno en vez de otra función periódica cualquiera es 
que la ecuación del movimiento se resuelve más fácilmente cuando 
se trabaja con una de esas dos funciones. Cuando esta nueva fuerza 
se tiene en cuenta, la segunda ley de Newton se convierte en 


mX"(t) = fuerza total ejercida sobre el cuerpo 


—sX(t) + F, cos pt (te RG). 


Generalmente esto se escribe de modo que la función desconocida X 
quede en el miembro de la izquierda: 


mX"(t) + sX(t) = F, cos pt (tER¿). 
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Ejercicio 1 
(3 minutos) 


Ejercicio 2 
(5 minutos) 


31.3 


31.3.1 


Discusión 
* * 


Definiciones 1 
* y 


(continúa en la página 28) 
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Solución 31.2.4.1 Solución 31.2.4.1 
En este caso, 

X(0)=x. y  X(0)=00, 
de modo que 


v 
X :tH—> x COS wt + sent (te R¿). |] 
w 


Solución 31.2.4.2 Solución 31.2.4.2 


El número de oscilaciones por segundo es 1/T, lo cual es igual a w/2r. 
Para aumentar ese número son apropiadas las alternativas (b) y (c) 
(aumentar s o disminuir m). 


Esto también se puede comprender fácilmente si se toma directamen- 
te una base intuitiva en la mecánica: una masa liviana, sometida a 
la acción de un muelle fuerte, se mueve rápidamente, mientras que 
una masa pesada adherida a un muelle débil se moverá lentamente. 


(viene de la página 27) 


Esta ecuación del movimiento, con una escogencia particular de p, es 
la base de nuestro estudio de la resonancia. Después de resolver el 
caso general como un problema matemático, interpretaremos su solu- 
ción en términos del comportamiento del sistema físico que representa. 


31.3.2 Solución de la ecuación del movimiento 31.3.2 
Frente al problema de la solución de la ecuación del movimiento Técnica 
* h 


mX"(t) + sX(t) = F, cos pt (teR¿), 


podemos, una vez más, recurrir a la lista de preguntas de Polya. 
COMPRENSION DEL PROBLEMA 


¿Cuál es la incógnita? Un conjunto de funciones reales de la 
forma f¿Ro —>R. 


¿Cuáles son los datos? m, s, Fo y p. 
¿Cuál es la condición? El conjunto debe contener todas las solu- 


ciones de la ecuación dada. 


TRAZADO DE UN PLAN 
¿Conoce algún problema relacionado con éste? 


Precisamente acabamos de trabajar en el problema relacionado con 
éste que consistía en la solución de la ecuación 


mX"+sX=0, 
Le. 
mX"(t) + sX(t) = 0. 


Aquí hay un problema semejante y que se ha resuelto antes. ¿Podría 
usted usarlo? 
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Sabemos que el conjunto solución de mX" + sX = O es el espacio 
vectorial generado por las dos funciones 
t—— COS 0wt, tt——senot. 


La ecuación mX" + sX = O es un caso especial (con Fy = 0) de la 
que estamos examinando ahora y, por tanto, nuestro problema consiste 
en generalizar la solución de mX” + sX = O al caso en que el miem- 
bro de la derecha es una función diferente de cero. 


¿Podría resolver una parte del problema ? 


Podríamos intentar encontrar solamente una parte del conjunto solu- 
ción en vez del conjunto completo; esto es, una solución particular de 
nuestra ecuación 


mX"(t) + sX(t) = F, cos pt. 


Sucede que esta ecuación tiene una sencilla solución particular. ¿Ve 
el modo de obtenerla? 


Ejercicio 1 Ejercicio 1 
> , , (3 minutos) 
Demostrar que, si p? 4 s/m, entonces existe un número k tal que 
X :tH——> k cos pt (teR¿) 


es una solución de 


mX"(t) + sX(t) = F,¿cos pt. ] 
. Ejercicio 2 Ejercicio 2 
A z (3 minutos) 
Encontrar una solución particular de 
X"(t) + X(t) =sen pt, (teR¿) 
supuesto que p?* A 1. E 
Volviendo a nuestro problema, tenemos una solución particular de la Técnica 
ecuación diferencial, pero es sólo una solución parcial de nuestro pro- *o* 
blema porque puede haber otras soluciones de la ecuación. Para en- 
contrarlas continuaremos con las preguntas de Polya. 
¿Podría usted reformular el problema? 
La ecuación se puede escribir de varias maneras. Por ejemplo, pode- 
mos usar la notación de función en vez de la notación de imagen. En- 
tonces, 
mX"” + sX = tH—>F,Cos pt 
O, más concisamente, 
mX"+sX=g 
donde g:t-—>Fo, cos pt. Una forma más concisa aún es 
LUX) = 8, 
donde L es el operador lineal que aplica Xa mX”" + sX. 
Ya antes hemos encontrado ecuaciones de esta última forma. En la 
Unidad 26, Algebra lineal III, Sección 26.2.2, estudiamos sistemas de 
ecuaciones lineales de la forma 
Ax =b, (continúa en la página 30) 
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Solución 1 


Llevando a la ecuación diferencial la forma sugerida de X obtenemos 
la condición 


m(— p?k cos pt) + sk cos pt = F, cos pt 
Le: 
(—mp?*k + sk) cos pt = F, cos pt, 
que se satisface si 
—mp*k + sk = Fo. 


Puesto que p? 4 s/m, esta ecuación se puede resolver para k, de 
manera que 


ko Fo 
s — mp? 


La solución particular pertinente de la ecuación es, por consiguiente, 


Fo 


Solución 2 


La ecuación dada es igual a la anterior, pero tiene un seno en lugar 
de un coseno. Esto sugiere que probemos 


X :tt—> ksenpt (te Ry) 


como una posible solución, y la correspondiente sustitución en la ecua- 
ción diferencial muestra que esta función verifica la ecuación, donde 


1 
k = —. 
1 — p? 
Una solución particular de la ecuación es, pues, 
1 
xi Jsenos (teR¿). A 


(viene de la página 29) 


donde x y b son vectores columnas y A es una matriz n X n. Encon- 
tramos que, si A es no singular, entonces Ax = b tiene solución 
única; pero que, si Á es singular, entonces Ax = b tiene muchas so- 
luciones, supuesto que tenga alguna. En el último caso, podemos en- 
contrar muchas soluciones así: si z es un elemento del núcleo de la 
aplicación representada por A (esto es, una solución de Ax = 0) y si 
xo es alguna solución cualquiera de Ax = b, entonces xo + z es 
otra solución de Ax = b; en efecto, tenemos que 


Á(Xo + 2) = Axo + Az=b+0=b. 


No es difícil demostrar que toda solución de la ecuación Ax = b es 
de la forma xo = z, donde z pertenece al núcleo. (Esto último lo de- 
mostraremos en el texto.) 


DESARROLLO DEL PLAN 
Ahora podemos desarrollar el tercero de los cuatro principales pasos 
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Solución 1 


Solución 2 


que señala Polya en el proceso de la solución de un problema. Sabe- 
mos que la solución general de la ecuación matricial 


AX'= b 
es 


x= (alguna solución particular de Ax = b) 
+ (el elemento general del núcleo de A). 


También sabemos que nuestra situación aquí es efectivamente la mis- 
ma: la función X pertenece a un conjunto de funciones que forman 
un espacio vectorial (correspondiente al espacio vectorial al cual per- 
tenece x), y el operador L es una trasformación lineal en este espacio 
«vectorial (correspondiente a la trasformación lineal representada por 
A). Por consiguiente, podemos esperar que encontraremos la solución 
general de la ecuación diferencial L(X) = g de una manera análoga: 


X= (alguna solución particular de L(X) = g) 
+ (el elemento general del núcleo de L), 


donde el núcleo de L está definido como el conjunto de funciones que 
están aplicadas a O por medio de L. Una expresión general de todos 
los elementos del núcleo de un operador diferencial L recibe, general- 
mente, el nombre de función complementaria de L porque “comple- 
menta” la solución particular para dar la solución general. 


Ejercicio 3 
Para el caso en que 


et —> PF ¡00s pt 


LiX=>mX" + 5X, 


hemos encontrado una solución particular de L(X) = g, y también 
* la solución general de L(X) = 0. Use esta información para escribir 
una fórmula para la solución general de la ecuación diferencial 


mA sx == >AP0 ¡COSipt, 
donde 


le 
Pn-. e 


UNA MIRADA ATRAS 


Habiendo descubierto un método que parece resolver el problema, pue- 
de suceder que usted sienta la tentación de no continuar en este te- 
ma y quiera pasar lo más rápidamente posible a otra cosa. Es un error 
dejarse arrastrar por esta tentación; es posible que usted pierda el 
esfuerzo empleado en encontrar una solución. La supuesta solución 
puede estar incompleta en alguna forma, o ser incorrecta; además, 
puede haber alguna lección útil que se pueda aprender de ella. Por 
esta razón, el cuarto y último paso de Polya para resolver problemas 
nos dice que examinemos la solución obtenida y que, en particular, 
comprobemos los resultados, Esto es, nos volvemos al problema origi- 
nal y verificamos que realmente tenemos la solución, nada más ni na- 
da menos. 
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Definición 1 


* $ 


Ejercicio 3 
(3 minutos) 


Técnica 
* $ 


MB 31.3.2, 31.3.3 


En el caso presente, nuestro problema consiste en encontrar la solu- 
ción general (es decir, el conjunto solución) de la ecuación diferencial: 
mX"(t) + sX(t) = F¿ cos pt (teR¿), 


y hemos llegado a la supuesta solución: 
F, 
X(t) = (E) cos pt + a.cos wt + Bsenot, 
s— mp 


donde a y $ son números reales arbitrarios. 


Los dos ejercicios siguientes, considerados conjuntamente, piden que 
se compruebe que esta es verdaderamente la solución general. En am- 
bos ejercicios debe usted suponer que p? X s/m. 


Ejercicio 4 Ejercicio 4 
S (3 minutos) 
Compruebe que toda función de la forma dada anteriormente es una 


solución de la ecuación diferencial, cualesquiera que sean los valores 


de a y 6. | 
Ejercicio 5 Ejercicio 5 
(5 minutos) 
Compruebe que toda solución real de la ecuación diferencial tiene la 
forma dada anteriormente. 
La sugerencia de Polya que quizá es más útil aquí es: “Aquí hay un 
problema relacionado con usted y que ha sido resuelto anteriormente. 
¿Puede usted usarlo? Mn 
31.3.3 Interpretación de la solución 31.3.3 
Nos falta interpretar la solución general que hemos encontrado y que. Tema principal 
corresponde a la ecuación pios 
s 
mX"(t) + sX(t) = F¿ cos pt (» E 2. 
m 
Nuestra solución tiene la forma 
Fo ) + 
X :tt > | ——=]| cos pt + acos wt + Bsenwt (teRg). 
s— mp 


Como ya hemos visto, el miembro de la derecha es la suma de 
(i) una solución particular: 


Fo 
tt |] cos pt teRó » 
: 2 sp (teRy) 
donde la imagen es proporcional a la fuerza aplicada F, cos pt , 
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y no contiene constantes arbitrarias (excepto las que describen 
la fuerza aplicada misma); 


3 
(ii) una función complementaria: 


tacos wt + Psenot (te RG), 
que no depende de la naturaleza de la fuerza aplicada y que con- 


tiene dos constantes arbitrarias como en el caso de las vibracio- 
nes libres. 


A causa de las dos constantes arbitrarias a y $, la expresión corres- 
pondiente a X (t) produce una familia completa de funciones que des- 
criben posibles movimientos del sistema, uno por cada par de números 
(a, 8). Para escoger el elemento particular de esta familia que corres- 
ponde a un movimiento particular del sistema, necesitamos dos datos 
informativos; generalmente éstos son la posición y la velocidad del 
cuerpo vibrante en el momento en que se inicia el movimiento. 


Ejemplo 1 

Supongamos que la partícula parte del reposo desde el origen; esto es, 
X (0) = 0 (sale del origen), 
X'(0) = 0 (parte del reposo). 


La solución general da 


F F 
xo) =| 2 5) 0050 + acos0 + pseno= [2] 4 a 
Ss — mp Ss — mp 


—pF, 
X"(0) = [222 seno — awsen0 + fwcosO0 = Bo, 


y usando las condiciones X (0) = X'(0) = 0, obtenemos 


Por consiguiente, la solución particular de la ecuación adecuada a la 
información dada es 


F Fl 
Xx H> SS cos pt — [0] cos Wwt (te R¿). 


La siguiente gráfica ilustra el caso en que Fo =s=m=1,p=. 
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Ejemplo 1 


(continúa en la página 35) 


Solución 31.3.2.3 
La solución particular de L(X) = g es 


Fo 


tt— 
le 


s) cos pt (teR¿). 


(Véase el ejercicio 1.) 


La función complementaria, es decir, el núcleo de L, es la solución 


general de mX” + sX = O, que ya encontramos en las Secciones 
31:22 y 3: 
tacos wt + Bsenwt (teR¿). 


Reuniéndolas encontramos la solución general de L(X) = g: 
F p 
— (E) cos pt + a.cos wt + Bsenwt (teR¿), 
s— mp 
supuesto que p? 4 s/m. L 


Solución 31.3.2.4 


La derivación de la supuesta solución da 


F, 
Xt) = — [PE sen pt — awsenwt + fw.cos wt, 
s— mp 
p?F 
X"(t) = — [e] cos pt — aw? cos wt — Bu?*senot, 
s— mp 
de donde 
—mp? F, cos pt + sF, cos pt 
mX"(t) E sX(t) sn | Prensas = Fo Cos pt. 
s— mp 
como se pedía. ll 


Solución 31.3.2.5 


Ya hemos notado la analogía entre la ecuación que consideramos 
aquí, que tiene la forma L(X) = g, y la ecuación matricial Ax = b. 
El problema “relacionado con usted” es el demostrar que, si xo es 
una solución particular de Ax = b, entonces la solución general de 
Ax =b es 


X = Xo + (el elemento general del núcleo de A) 


La demostración se dio en la Sección 26.2.2 de la Unidad 26, Algebra 
lineal III, y la repetiremos aquí. Sea xo la solución particular de 
Ax = b, y sea x otra solución cualquiera, de manera que 


AX = b 


ÁXo = b. 


Sustrayendo y, además, por ser A una trasformación lineal, tenemos 
que 


Á(X — Xo) =0, 


y, por tanto, x — xo pertenece al núcleo de A, esto es, x es de la 
forma 


xo + un elemento del núcleo de A 
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Solución 31.3.2.3 


Solución 31.3.2.4 


Solución 31.3.2.5 


Así, cualquier solución de Ax = b es de esta forma; en otras pala- 
bras, la solución general tiene la forma 


xo + (el elemento general del núcleo de A) 


Para demostrar aquí el resultado correspondiente, necesitamos algo 
más que un cambio de notación. La ecuación diferencial dada se pue- 
de escribir como 


UX) = g 
donde 
L:X—>mX"+sX y  g:ttH——>F,cospt. 


Ya hemos visto que una solución de esta ecuación es la función 
Fo 
tt | ——5 | cos pt teRi 
: a ma) Pp ( 0 iS 


que denotamos por X/. Entonces, si X es cualquier solución de la 
ecuación diferencial, tenemos que 


LX) =8 


LUX) =g2. 


Puesto que L es una trasformación lineal, la sustracción nos da 
L(X — Xo) = O, de manera que X — Xo debe pertenecer al núcleo 
de L, esto es, al conjunto solución de la ecuación L(X) = O. Ya he- 
mos visto en la Sección 31.2.3 que este conjunto comprende todas las 
funciones de la forma 


tacos wt + fBsenot  (teR¿); 
por consiguiente, X— Xoy debe ser de esta forma y, por tanto, tene- 
mos, por definición de Xo, 


F, 
X:it— (E) cos pt + acos wt + fsenot  (teR¿), 


para algún «a y algún $. Así, toda solución X de L(X) = g tiene 
la forma establecida. a 
(viene de la página 33) 

Ejercicio 1 


Si la partícula sale de una posición xy con una velocidad vo, ¿cuál 
es la fórmula de X (t) que describe el movimiento correspondiente? MW 


Ejercicio 2 
Si la ecuación del movimiento es 
X"(t) + X(t) =senjt  (teRj), 


y la partícula sale del reposo desde el origen, encuentre la fórmula de 
la elongación X (t). 


(Véase el ejercicio 31.3.2.2.) L 


Las gráficas del ejemplo 1 y de la solución 2 muestran que el movi- 
miento en las vibraciones forzadas es, probablemente, más complica- 
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Ejercicio 1 
(4 minutos) 


Ejercicio 2 
(5 minutos) 


Discusión 
* * 


do que en las vibraciones libres. Esto no es más que lo que debemos 
esperar porque la ecuación diferencial correspondiente a las vibracio- 
nes forzadas es mucho más complicada. La gráfica pierde el carácter 
de curva sinusoidal de período 21/w que tiene la de las vibraciones 
libres; en general, consiste en una superposición de esta onda sinu- 
soidal y otra, cuyo período está determinado por la función que da la 
fuerza aplicada. La superposición de estos dos movimientos de perío- 
dos diferentes nos lleva a movimientos bastante complicados (que no 
tienen que ser necesariamente periódicos). 


Ejercicio 3 
(i) ¿Cuál es el período de la función 
tH——> COS pt (te Ry? 
(ii) ¿Cuál es el período de la función 
X:tt—¿senit — ¿sent  (teR¿) 


que encontramos en la solución 2? | 


Como indicamos al comenzar este texto, una importante aplicación de 
la matemática de las vibraciones la tenemos en el fenómeno de la re- 
sonancia. Por resonancia entendemos la tendencia de las oscilaciones 
de un sistema vibrante a alcanzar amplitudes muy grandes cuando 
la fuerza que actúa sobre él oscila con un período que es igual al perío- 
do de sus vibraciones libres o que está muy cerca de él. Nuestro mo- 
delo explica cuantitativamente este fenómeno. Supongamos, para fijar 
ideas, que la fuerza aplicada es F, cos pt y que sale del reposo desde 
el origen en el tiempo 0. Entonces, por la solución 1, el movimiento co- 
rrespondiente viene dado por 


F 
x:0— Eos pt — cos wt) (te RG). 
Ss — mp 


Los períodos de las vibraciones libres y de la fuerza aplicada son 2r/w 
y 27/p. Si esos períodos son casi iguales, tendremos que p está cerca 
de w lo cual implica que mp? está cerca de s y, por consiguiente, el 
denominador, (s — mp?), es pequeño. En consecuencia, el factor 


F , ; 
|, E) es grande y, por tanto, X(t) será grande, excepto en esos 
Ss — mp 


momentos (como en el tiempo inicial) en los cuales los dos cosenos se 
cancelan o casi se cancelan. Los diagramas muestran las gráficas de 
X, con Fo = m = s = 1, para dos valores diferentes de p. Cuanto 
más cerca esté p de w = 1, tanto mayor es la amplitud máxima. 


pul 


jo] 
tl 
aia 
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Ejercicio 3 
(2 minutos) 


Tema principal 


ramplitud; 


Para obtener una estimación sencilla de la máxima amplitud, supon- 
gamos que ésta ocurre en un valor de t tal que uno de los dos términos 
cos wt y cos pt está cerca de +1 y el otro está cerca de — 1. Enton- 
ces, |cospt — coswt|= 2 y, por tanto, la amplitud es, aproximada- 
2Fo 
s— mp 


mente, Esta cantidad se representa gráficamente, con p 


34" 


como variable independiente, en este diagrama: 


Vemos que la máxima amplitud estimada puede, en realidad, ser muy 
grande cuando p está cerca de «u. 


¿Qué sucede si p llega a ser, realmente, igual a w, de manera que la 
fuerza aplicada tenga exactamente el mismo período que las vibracio- 
nes libres? La solución que hemos venido usando se basa en el su- 
puesto que p? + s/m, de manera que no comprende este caso. La 
ecuación del movimiento es, ahora, (con Fy = s = m para mayor 
sencillez) 


X"(t) + X(t) = cost 


y no es posible lograr que se verifique la ecuación si tomamos X(t) = 
k cost, donde k es un número, porque el miembro de la izquierda de 
la ecuación nunca será, en ese caso, igual a cost sino que se anulará 
para todos los valores de k. 


Hay varios métodos para encontrar soluciones particulares en estos 
casos excepcionales. Uno de ellos se indica en el apéndice. Sin em- 
bargo, todo lo que nos interesa aquí es el resultado, que es, precisa- 
mente, que la función 
Xx: —Hsent (teR¿) 

es una solución particular de la ecuación anterior. En el diagrama 
mostramos la gráfica de esta función. Nótese que comienza de mane- 
ra muy semejante a la gráfica en que p = 8/9. La diferencia con- 
siste en que ésta continúa creciendo indefinidamente: nunca alcan- 
za un máximo. 
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(continúa en la página 39) 


Solución 1 


Comenzando como en el ejemplo, tenemos ahora que 


F 
[7] +2=x y  Bw=00, 


de modo que 
F v 
A = Xy — [hal y B= pen 


s — mp? 


El movimiento corresponde a la función 


Fo Vo 
———3|/008 0 + =senot 
mp w 


Xu— pe cos pt + |xy— 
s — mp? s— 


(te R¿) n 


Solución 2 


En el ejercicio 31.3.2.2 encontramos una solución particular de la 
ecuación 


X"(t) + X(t) =senpt, 


a saber 
1 4 
>El 7 7|sen pt (te Ro). 
=p 
Aquí tenemos la misma ecuación con p = 3, de manera que una so- 


lución particular es 
X :tt—Ísenjt (teR¿). 


La solución general de la ecuación del movimiento, obtenida de la 
adición de esta solución particular al elemento general del núcleo 
del operador 


e o E O A 
es 


X:tt——sentt + acost + Bsent — (teRj¿). 


Para usar las condiciones iniciales observamos que 


X(0) = Ísen0 + a.cos0 + fsen0 = a 


X"(0) =3c0s0 — asen0 + Bcos0 =34 + Pf, 


de manera que a = 0 y 8 = — 3, si la partícula sale del reposo desde 


el origen. Sustituyendo estos valores de a« y f en la solución gene- 
ral anterior obtenemos la solución buscada 


X :tt——>$senjt — ¿sent (te R¿). 
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Solución 1 


Solución 2 
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Solución 3 . Solución 3 


(1) 


tr» Cos pt 


Vemos que la región de la gráfica comprendida en el intervalo 
2x1 ; o 

[0.2] se repite una y otra vez, y que ningún intervalo menor 

que ese tiene esta propiedad. El período es, por consiguiente, 27 /p. 


(ii) En la solución 2 no trazamos suficientemente la gráfica y, por 
consiguiente, no se puede observar el período pero esto se puede 
suplir fácilmente. 


El período del término sen 3t es 4r, y el de sen t es 2r. Por tan- 
to, el período de la función completa es 4r. Es 


(viene de la página 37) 


de 
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Ejercicio 4 
Verifique que 
X :t—Htsent 


es una solución particular de 
X"(t) + X(t) = cost. E 
Ejercicio 5 
Escriba la solución general de 
X"(t) + X(t) = cost, 
y escoja la solución particular que describe el movimiento de una 


partícula que, en el tiempo 0, tiene una velocidad vo y una elonga- 
ción xo. a 


Ejercicio 6 


Es obvio que no es realmente posible que las vibraciones aumenten 
la amplitud arbitrariamente. La predicción de nuestro modelo mate- 
mático es que se podrá en el caso en que p = 1. ¿Qué característica 
de la situación física, omitida en nuestro modelo matemático, debe 
incluirse si queremos que el modelo dé una descripción más realís- 
tica del caso en que p = 1? | 


En realidad, la ecuación del movimiento que hemos dado en el texto 
es una buena aproximación solamente cuando se aplica a pequeñas 
elongaciones. Si las elongaciones son grandes, tendremos que usar 
modelos matemáticos más complicados. 


No obstante, el modelo que hemos estudiado predice correctamente y 
describe el crecimiento inicial de las oscilaciones, antes de que co- 
miencen a ser muy violentas. 


Posdata 


En la vida real hay algunos ejemplos sorprendentes de resonancia. 
Algunos son útiles y otros son desastrosos. Un ejemplo de los prime- 
ros es un circuito eléctrico sintonizado para recibir señales de un tipo 
particular como las de un trasmisor de radio. El circuito (el radio- 
rreceptor) responde únicamente al estímulo de una frecuencia parti- 
cular: la frecuencia de resonancia. 


Por otra parte, en 1940 el puente de suspensión que estaba sobre el 
desfiladero de Tacoma en Norteamérica fue destruido por el viento 
que lo hizo vibrar a su frecuencia natural. Hay una película tomada 
por un observador que muestra las oscilaciones desde que comenza- 
ron hasta que provocaron la caída del puente. Algunos instrumentos 
musicales ilustran la resonancia; si cerca de un violín se toca una 
nota adecuada, la cuerda correspondiente comenzará a vibrar por 
simpatía. Ligeramente diferente, pero de la misma clase, es el caso 
que se cuenta de la famosa soprano Nellie Melba que era capaz de 
romper una copa de cristal cantando una nota suficientemente alta. 


Otro ejemplo se ve en los soldados que “rompen el paso” cuando 
marchan sobre un puente colgante. Si no lo hicieran, la fuerza ejer- 
cida sobre el puente podría coincidir con la frecuencia natural, y la 
resonancia podría destruir el puente. 
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Ejercicio 4 
(3 minutos) 


Ejercicio 5 
(4 minutos) 


Ejercicio 6 
(5 minutos) 


Discusión 
* * 


Canto de amor de dos muelles 


“Amo el modo en que te enroscas”, 
dijo la bobina al 
muelle. El muelle 
dijo: “Las oscilaciones 
son ahora la última 
moda”. “Hay un 
período de tiempo, que 
con frecuencia ocurre en primavera, 
cuando las fuerzas externas 
resuenan y te arrastran hasta 
el bamboleo.” “Pero, en 
otoño”, dijo entonces 
con esquivez, “cuando la 
fricción está 
creciendo, la 
lluvia y 
la humedad 
hacen que 
la amplitud 
decrezca”. 
“Así que por ahora 
viviremos en 
armonía, nuestras 
funciones 
complementando, 
y cuando las osci- 
laciones mueran, 
podremos 
entonces 
descansar 
en 


31.4 APENDICE (no forma parte del curso) 
Una solución particular de X"'(t) + X(t) = cos t 


Para delimitar las ideas, busquemos una solución particular que co- 
rresponde a un movimiento que parte del reposo desde el origen. Usa- 
remos un proceso de límite en el cual p se acerca a 1, aplicado a la 
correspondiente solución de la ecuación 


X"(t) + X(t) = cos pt (te R¿). 
Según se mostró en el texto, con 


o 5=m=1 
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(continúa en la página 42) 


Solución 4 
Si 
X(t) = Htsent, 
entonces 
X'(t) = ¿sent + Htcos t, 
X"(t) = 3 cost + ¿cost — ¿tsent, 
y, por tanto, 


X"(t) + X(t) = cost como se pedía |] 


Solución 5 
La solución general es 


X :tt——dtsent + acost + fsent, 


de donde 

X(0) = a 

AO) 8: 
Si la partícula tiene en el momento t = 0 una velocidad vo y una 
elongación xo, tenemos que a = xo, $ =U0, y 

X :tt——htsent + xpcos t + vosent. 7 
Solución 6 


Las principales idealizaciones que hemos hecho en el modelo son: 


(i) no se ha tenido en cuenta la fricción; 
(ii) se ha supuesto que, cualquiera que sea la elongación, la fuerza 
sobre la partícula es —sX(t) + (fuerza ejercida); 
(iii) las oscilaciones no llegan a alcanzar el soporte. 


Así, pues, suponemos que la partícula no golpea ningún otro objeto 
(por ejemplo, su soporte), que las propiedades del muelle no cambian 
con el tiempo (en particular, no se rompe) y que la fuerza que hay en 
el muelle varía linealmente con la elongación X(t) sin que importe 
cuán violentamente se estire o comprima. 15 


(viene de la página 41) 


la solución de esta ecuación con 


X(0) = X(0)=0 
es 
cos pt — cost 
Ly AA (te R¿). 
dB 

Esto sugiere que 

X :t—> X(t) = lím [a (teR¿) 

p-=1 1 — Y) 


puede ser una solución de X”(t) + X(t) = cost. Para comprobar es- 


42 


MB 31.4 


Solución 4 


Solución 5 


Solución 6 


ta conjetura, evaluamos el límite. Esto se puede hacer mediante algu- 


nas manipulaciones basadas en la identidad* 


cos (u — v) — cos (u + v) = 2senusenv 


la cual, con u y v escogidos de modo que u+u=tyu-=—u 


nos da 


l+ 1 — 
cos pt = cost = 2sen]| esco ul 


Encontramos que 


lím cos pt — cost 
p>1 1 — DE 
2sen de tsen a 
= lím 2 * 
p>1 (1 + pd —p) 
1 1-— 
2sen] al sen > al 
= lím x lím 
p=1 1 + p>1 fs A 


(por la regla del producto de límites (véase la Unidad 7)). 


= DE, 


2sent t , senx , 
= x1x>3 puesto que lím—— = 1, según 
2 2 EA IO 
tsent se ha demostrado en el apéndice II 
ar, Y de la Unidad 12, Diferenciación 1 


Esta es la forma que se le pedía en el texto que considerara para X (t). 


*Este no es el método más sencillo pero, sin embargo, es el que requiere menos co- 


nocimientos previos. 
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M 100 — UNIDADES DEL CURSO BASICO DE MATEMATICAS 


0J3050NAO0NA 


Funciones 

Errores y exactitud 
Operaciones y morfismos 
Diferencias finitas 

NO HAY TEXTO 
Desiguales 

Sucesiones y límites 1 
Computación Il 

Integración I 

NO HAY TEXTO 

Lógica I — Algebra de Boole 
Diferenciación 1 
Integración II 

Sucesiones y límites II 
Diferenciación II 
Probabilidad y estadística 1 
Lógica Il — Prueba 
Probabilidad y estadística II 
Relaciones 

Computación IH 
Probabilidad y estadística HI 
Algebra lineal 1 

Algebra lineal II 
Ecuaciones diferenciales 1 
NO HAY TEXTO 

Algebra lineal II 

Números complejos 1 
Algebra lineal IV 

Números complejos Il 
Grupos 1 

Ecuaciones diferenciales II 
NO HAY TEXTO 

Grupos II 

Sistemas numéricos 
Topología 

Estructuras matemáticas 
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